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We articulate the fact that the loop quantum gravity description of the quantum states of black 
hole horizons, modeled as Quantum Isolated Horizons (QIHs), is completely characterized in terms of 
two independent integer- valued quantum 'hairs', viz,, the coupling constant of the quantum SU(2) 
Chern Simons theory describing QIH dynamics, and the number of punctures produced by the bulk 
spin network edges piercing the isolated horizon (which act as pointlike sources for the Chern Simons 
fields). We demonstrate that the microcanonical entropy of macroscopic (both parameters assuming 
^ ■ very large values) QIHs can be obtained directly from the microstates of this Chern-Simons theory, 

' using standard statistical mechanical methods, without having to additionally postulate the horizon 

as an ideal gas of punctures, or incorporate any additional classical or semi-classical input from 
general relativity vis-a-vis the functional dependence of the IH mass on its area, or indeed, without 
having to restrict to any special class of spins. Requiring the validity of the Bekenstein-Hawking area 
law relates these two parameters (as an equilibrium 'equation of state') and consequently imposes 
restrictions on the allowed values of the Barbero-Immirzi parameter. The logarithmic correction to 
the area law obtained a decade ago by R. Kaul and one of us (P.M.), ensues straightforwardly, with 
precisely the coefficient —3/2, making it a signature of the loop quantum gravity approach to black 
'^j ' ■ hole entropy. 
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^ . I. INTRODUCTION 

m 
in 

Perhaps the clearest description of the quantum states of generic (extremal or non-extremal) four dimen- 
sional black hole horizons in the presence of matter and/or radiation (which do not cross the horizons) is 
in terms of the quantization of the classical Isolated Horizon (IH) phase space which yields the kinematical 
£T) • Hilbert space for the Quantum Isolated Horizon (QIH) [1, 2] derived from loop quantum gravity (LQG). In 
this description, the physical states of a QIH are the gauge singlet states of an SU(2) Chern Simons (CS) 
theory coupled to the punctures endowed with SU(2) spins deposited on the QIH by the intersecting edges 
of the bulk spin network describing the bulk quantum geometry. In other words, the QIH is taken to have 
^ , the topology of a two-sphere with punctures carrying spins induced by the floating lattice known as the spin 
network whose edges carry spins. Such configurations are expected to arise as solutions of the full quantum 
dynamical equations describing bulk quantum geometry, in particular the quantum Hamiltonian constraint. 
Whether or not this actually happens remains a question for the future. In this situation, it is worthwhile 
to investigate properties of the kinematical Hilbert space to see how the issue of black hole entropy can be 
addressed. 

There are a plethora of approaches to black hole entropy, many semiclassical ones with diverse claims for 
the black hole entropy having logarithmic corrections to the area law [3, 6-10], with various coefficients - 
some postive and some negative. In all these computations, one usually computes the entropy of quantum 
'matter' fields (including gravitons) coupled to the classical background spacetime metric of a black hole. 
Even in approaches where the classical metric is 'integrated over', the functional integral over metrics is 
invariably saturated by the classical black hole metric. In other words, non-perturbatively large quantum 
fluctuations of the spacetime - which cannot meaningfully be separated into a classical background metric 
and its fluctuations - are usually ignored in these computations. Thus, these computations do not take 
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into account the entirety of quantum states corresponding to quantum spacetime fluctuations of a black 
hole spacetime. Rather, they account for non-gravitational states entangled with the horizon, as discussed 
in ref. [11, 12]. In contrast, the LQG approach focuses on the quantum states describing the quantum 
geometry of the horizon, without assuming any entanglement with quantum matter states in its vicinity. 
The classical metric of the black hole plays no role in this approach. The LQG computation thus yields what 
one may call the gravitational (or 'spacetime') entropy of a black hole, as opposed to the non-gravitational 
entanglement entropy computed in other approaches, which depends rather directly on the classical black 
hole metric in all cases. Indeed, the two together would most likely constitute the total entropy of a quantal 
black hole spacetime, so that the results of both directions of computation are complimentary rather than 
competitively comparable. In other words, the entanglement entropy considered in the other approaches 
must be an additional contribution to the horizon entropy, over and above the spacetime entropy computed 
within the LQG approach. 

The description of a QIH in terms of CS states coupled to bulk quantum geometry is itself a radical departure 
from a the general relativistic description of a classical event horizon. Of course, a QIH does not radiate 
or accrete matter/radiation, and can only be a part of a more complete description of a radiant black hole 
comprising of Trapping [18] or Dynamical [14] horizons. Recently, such horizons have been argued to emit 
Hawking radiation [15] with the standard black body spectrum, raising hopes that a deeper understanding 
of black holes within the LQG approach might be around the corner. However, there are semi-classical 
assumptions that appear to be necessary to supplement the premises of LQG, which manifest a lingering 
incompleteness in basic understanding of the quantum geometry of IHs. 

The object of the present paper is to articulate a fact that has already been implicit in recent literature [19], 
however, without the extra semi-classically motivated assumptions used earlier : that a complete description 
of QIH is possible with two integer- valued quantum parameters ('hairs'), namely the CS coupling constant k 
related to the classical horizon area, and the number N of punctures to which the CS fields couple. Indeed, 
the microcanonical entropy can be obtained directly from the formula derived fourteen years ago by one of 
us (P. M.) together with R. K. Kaul [20] for the total number of SU(2) singlet states coupled to a fixed 
set of spins j\, j'2, Jn , ji & [1/2, fc/2]Vi = 1,...N. Summing this formula over all spins {ji} and using 
the Multinomial expansion of elementary algebra, this degeneracy can be expressed in terms of sums over 
spin configurations, i.e., number of punctures rij € [0, N] for spins j = 1/2, k/2. Using standard tenets of 
equilibrium statistical mechanics, one looks for the 'Most Probable Spin Configuration' by maximizing the 
microcanonical entropy subject to the constraints [21] of a fixed N and fixed large QIH area differing from 
the classical area only to 0(£p), where £ p is the Planck length. Requiring that the resulting formula for the 
microcanonical entropy of a macroscopic (k ^> 1, N 1) QIH reproduce the Bekenstein-Hawking area law 
in terms of classical IH area A c i , imposes a restriction on the Barbero-Immirzi parameter 7 in the definition 
of k = A c i/4:TT'y£p, that this must lie within a specific interval on the real line [22]. Sub-leading corrections to 
the area law ensue naturally from the derived formula, especially the leading logarithmic (in classical area) 
correction with coefficient —3/2 found longer than a decade ago by Kaul and Majumdar [23] for a dominant 
class of spins, as we shall see in the sequel. These results have also been rederived recently in ref. [24-27] 

The important value-addition here is the direct link to the LQG formulation of QIHs in terms of SU(2) 
CS states coupled to punctures carrying spin, without any trappings of a semi-classical nature depending 
on classical metrical properties, nor any restriction to a specific class of spin values at the punctures. While 
many of the extant papers in the recent literature on QIH entropy within the LQG approach simply count 
states of an ideal gas of spins, without much allusion to the underlying theory of QIHs in terms of quantum 
CS states, our approach here has been to underline the CS theory underpinning, and to make the argument 
as self-contained as possible within the original LQG approach pioneered by Ashtekar and coworkers. This 
paper has some overlap with the recent review of Kaul [28] as far as some parts of the calculations are 
concerned, but the emphasis on direct link-up with the CS theory perspective and some of the details are 
different and are, hopefully, of inherent merit. 



II. PHYSICAL CHERN-SIMONS STATES OF THE QIH 

The computation of the microcanonical entropy formula of Kaul and Majumdar [20] proceeds from the 
expression for the number of conformal blocks of SU(2)k WZW model on a 2-sphere with marked points 
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(punctures) carrying spin which, in the seminal work of Wittcn [29] has been shown to give the dimensionality 
of the singlet part of the Hilbert space of SU(2) CS states onR®^ 2 coupled to punctures on the S 2 . Using 
this remarkable connection and the fusion algebra and the Verlinde formula, the degeneracy of the microstates 
is expressed as [20] 



fe+l gin Q7r(2ji+1) gin a-n(2j N + l) 
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To obtain the total number of conformal blocks, this expression must be summed over all possible spin 
values at each puncture. Before doing the sum over each spin, let us write down the quantum mechanical 
wave function for a QIH as a linear combination of the basis states represented by the spin network basis. 
Thus, a generic normalized wave function for the QIH, in this spin network basis, can be written as 

l*>= J2 5 0V- Jn) \h,--- Jn) (2) 

where (^1^) = 1 , ji € [1/2, fc/2]Vi = 1, ...,N. The latter requirement merely reiterates the fact that 
one must restrict to unitary representations of the SU(2)k current algebra. The total number of punctures 
N and the Chern-Simons level k being the two independent parameters of the theory [22], the number of 
microstates of a QIH characterized by given values of N and k can be written as 

ft(N,k)= ]T p(ji,--- ,3s) n(ji,--. ,j N ) (3) 

31,- " >3N 

where p{ji,--- = |c(ji,-- - , Jjv)| is the probability of finding the QIH in a particular eigen- 

state configured by the spin sequence (ji,-- - , X/v). It follows from the normalization condition that 
T,j 1 ,...,j N P{h, ■ ■ ■ Jn) = 1- 

Now, since we will apply the method of most probable distribution to find the microcanonical entropy of 
the QIH, it is convenient to switch over to the configuration basis. At first we write the wave function of the 
QIH in the configuration basis as = X){ s } c [{ s i}]l { s j})- Then the number of microstates for a QIH with 
a given N and k will be obtained by summing over all possible configurations(eigcnstates) and is written as 

n(iv ! fc) = ^cj[{ Sj }]o[{ Sj }] (4) 

{so} 

where J2{ S] } w [{ s il] = 2{ s } l c [{ s j}]| 2 = = li w [{ s il] being the probability to find the QIH in the 

particular eigenstate designated by the spin configuration {sj} and 

AM 



n[{ Sj }}= fj — -g[{s 3 }] (5) 



where j runs from 1/2 to k/2 as usual and 

The combinatorial factor in eq.(5) reflects the statistical distinguishablity of the punctures, a property 
inherited by the punctures in the quantization procedure of the classical IH due to the nontrivial holonomics 
of the Chern-Simons connection on the IH along disjoint closed loops about the punctures [2]. 



A. Multinomial Expansion : The Link 

It is worth mentioning that, even though eq.(5) resembles the formula for the microstates of an ideal gas 
of spins obeying Maxwcll-Boltzmann statistics, it has been derived directly from the formula (1) within the 
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present scenario of QIH. One should keep in mind that eq.(5) is just another form of eq.(3) being written 
in a different basis for convenience of the statistical formulation. Hence, for the sake of clarity, let us derive 
eq.(5) directly from eq.(3). 

First of all, since the punctures do not physically interact with each other, we can write p(ji, ■ ■ ■ , Jjv) = 
p{ji) ' ' ' pUn)- Before proceeding further, let us note this crucial point in slightly more detail. Although there 
is no physical interaction among the punctures (i.e. the spin at a puncture can take any value irrespective 
of spin values at other punctures) only those configurations giving rise to SU(2) singlet states are physical, 
the others being unphysical 'gauge' states eliminated by the quantum Gauss law constraint. The number of 
singlet states for a given ordered sequence of spins (jfi, • ■ • , j'at) can be written as a linear combination of 
delta functions [20]. Now, the formula for the microstates given by eq.(l) is directly derived from this delta 
function formula using Fourier transform. On the other hand, one can treat this delta function formula as a 
constraint which excludes the counting the extra 'gauge' degrees of freedom. Thus, in some sense, the spins 
at different punctures must be correlated to give rise to the physical states of the QIH. In fact this is the 
reason why the second order corrections to the entropy i.e. the logarithmic corrections [23] are not additive, 
which has been correctly pointed out in [30]. Hence, considering p(ji, ■ ■ ■ ,.]n) = p(ji) 1 1 ■ p(jn) is completely 
justified since the delta function formula being a constraint itself takes care of this virtual correlation among 
the punctures. 



Hence, using eq.(l) one can write eq.(3) in the following explicit form 
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Using Multinomial expansion, the above expression can be recast into the following form 
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(7) 



where ^[{sj}] = {p(j)} Sj ■ This is exactly eq.(4), which is nothing but eq.(3) written in the configuration 
basis. 



III. QUANTUM HAIRS AND MICROCANONICAL ENSEMBLE 



The next task is to define the microcanonical ensemble of QIHs by identifying and fixing the parameters 
of the theory. Right at the beginning, it has already been mentioned that k and N are the two independent 
parameters, whose values, once specified, determine the number of microstates of a QIH given by the eq.(7). 
Thus k and N completely characterize a QIH. To obtain the microcanonical entropy, one must maximize 
its expression with respect to the spin configurations, to determine the most probable spin configuration, 
subject to the restriction that the QIH area must equal the classical area upto 0(lp) and the number TV of 
punctures is fixed. 

The area eigenvalue equation for a particular eigenstate of the QIH in the configuration basis can be written 

as 

A\ { Sj }) = Svry^ s iVi(J+l) I {sj}) (8) 

3 
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Hence, the expectation value of the area operator for the QIH is given by 



(i) = (tt|4|tt) = 87T 7 4£ W [{ S ^]5>V7(J+l) = A cl ±0(f p ) (9) 

M j 

where A c i is the area of the classical IH closely represented by the QIH. Scaling the equation by 8irj£p and 
using (A)/8ir~f£l w A c i/^l 2 p = k/2 [2, 16], we obtain 

^c[{ Sj }]^ Sj V7(7 + T) = ^ (10) 

In this process we can also avoid the involvement of the ambiguous parameter 7 which will be ultimately 
fixed at the end by the usual argument of the validity of the BHAL. 

Apart from this, the expectation value of the number of punctures for the QIH is given by 

(AO = <*|JV|*) = = N (11) 

{sj} 3 

where TV can be considered as the operator for the number of punctures for a QIH. One should note that 
unlike the case of area the expectation value of N is exactly equal to the total number of punctures N. 
Now, since we are doing equilibrium statistical mechanics (thus stability is assumed) of QIH with large area 
(A c i ^> £p) and large number of punctures (N 1), for all practical purposes, we can neglect the fluctuations 
and also consider that the dominant contribution to the entropy comes from a most probable configuration, 
{s*} i.e. w[{s*}] ~ 1. Hence, the spin configuration {s*} must obey the following constraints 

d:J2 s *j = N ( 12a ) 

3 

C 2 :^4V7(7 + T) = | (12b) 

3 

Hence, we define a microcanonical ensemble of QIHs by assigning fixed values of k and N respectively. The 
obvious following step is the computation of the microcanonical entropy of a QIH characterized by fixed 
values of k and N. 



IV. MICROCANONICAL ENTROPY 



Having defined the microcanonical ensemble appropriately, we shall now derive the microcanonical entropy 
of a QIH. The microcanonical entropy of a QIH for given values of k and N is written as 



S MC = logn(N,k)~logn[{s* i }] 



(13) 



where we have set the Maxwell-Boltzmann constant to unity. Variation of log f2[{sj}] with respect to Sj, 
subject to the constraints C\ and C2, yields the distribution function for the most probable configuration 
{s*} which maximizes the entropy of the QIH. In other words, s* satisfies the variational equation written 
as 



6 logn[{ Sj }) -a^6 Sj -Xj2 hi VJti + V = 



(14) 



where i5 represents variation with respect to Sj , a and A are the Lagrange multipliers for C\ and C2 respectively. 
This yields the most probable distribution given by 



N exp 



OSj 



(15) 
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To proceed further we calculate (/[{s^}] explicitly using saddle point approximation in the limit k, N — > oo, 
which is appropriate for large black holes. First of all, we rewrite eq.(6) replacing the summation over o by 
integration as 
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where we have applied a change in the integration variable as air/(k + 2) = 9 and for which the limits follow 
with e = Tr/(k + 2). Now, for k — > oo, e — > 0. Hence, we can safely write 
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= - f cxp [G(6, k)} d6-- [ cxp [ln(cos 29) + G(9, k)} d9 (17) 
71 Jo n Jo 

where G(9, k) = J2j Sj log { b ' n ^ J n ^ 9 |. The above two integrations can be performed by the saddle point 
method. To begin with, it is straightforward to show that 
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which implies that the saddle point 9 ~ ( ' denotes partial derivative with respect to 9). Now, Taylor 
expanding G(9, k) about the saddle point 9q up to second order and applying the saddle point conditions 
from eqs.(18), we have 
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where we have used £ = 9 — 9q and the following limits 
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Evaluation of the integral yields 
1 
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The quantity a results from the second order approximation. Hence, while calculating a we can only use 
the results up to first order i.e. we use s* in place of Sj whose expression will be given by 



3* ~ N(2j + 1) cxp[-Av%- + l) - a] 



(21) 



which follows from the fact that g[{sj}] — ^ ]^[ - (2j + l) Sj neglecting the second order corrections, C being 
some constant. Now, using eq.(21) in eq.(12a) and eq.(12b), one obtains 

exp[a] = ^(2i + l)cxp[-AV7(7TT)] 



k/2 = N J2 VJU + l)(2j + 1) exp[-AVi(i + 1)-<t] 



(22a) 
(22b) 
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Using eq.(22a) and eq.(22b) one can show that k = ~2N (da / dX) and consequently a ~ %N (d 2 a / dX 2 ) . Now, 



^V(4 + «)/2 



Erf 



the quantity a — > oo for A^ — > oo. Hence, we can approximately write Erf 

Plotting the function Erf ny/a/2 with a one can see that the function attains a constant value for large 



a. Hence, we can take lim Q 
from eq.(20) it follows that 
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K, some constant. Therefore, 




Therefore, we have 



Ssj 



logg[{ Sj }\ = log(2j + 1) - + 1) 



(23) 



Hence, considering variation of a resulting from the inclusion of the quadratic fluctuations, the distribution 
for the most probable configuration given by eq.(21) gets modified into 



N(2j + l)exp[- A y/j (J+T) 
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(24) 



Now, using the most probable distribution given by eq.(21) or eq.(24) (result will differ by a constant only) 
we calculate the microcanonical entropy of a QIH for given values of fc and N using Stirling approximation 
and the result comes out to be 



Afc 

i>Aic = — 



Na - | log N - ^ log(d 2 cr/(iA 2 ) + 



(25) 



where A and a are explicit functions ofk/N which can be obtained by solving equations (22a) and (22b) [22]. 
Now, for a QIH (i.e. at equilibrium), A c i and TV are related by the equation A c i = S-K^i 2 N{da/d\) and 
also we have fc = A c i/4irj£p. Hence, the microcanonical entropy can be expressed completely in terms of 
the classical area of the IH or the mean area of the QIH as 
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where /(A) = A— (cr/4f )■ To retrieve the Bekenstein-Hawking area law(BHAL), one should fit 7 = f(X)/2n. 
One should note that, in this process, the value of 7 becomes known in terms of fc and N as A is completely 
solvable in terms of those parameters from the equations (22a) and (22b). This issue has been explicitly 
discussed and shown in [22] , although there is a subtle difference in the fitting of 7 as far as the present work 
is concerned. 



Now we calculate a explicitly in the limit fc — > 00 to study the nature of the function /(A). We have 

poo 

lim cxp[er] ~ / (2j + l)cxp -XyJ j{j + 1) dj - 1 



which yields a = log 1 



1) and da/dX = 4/[A(A 2 



/(A) = A 



2)]. Hence, it follows after a few steps of algebra 
A 2 



A ) log 
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(27) 



(28) 



Now, from information theoretic arguments it can be explained that a < [22], which in turn implies that 
A > 1. On the other hand, for /(A) to be real (for positive entropy), we must have A 2 < 2. Thus the allowed 
(positive) values of the parameter A get restricted in the range 1 < A < \[2. Also, one can check that the 
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corresponding range of the function /(A) is also given by 1 < /(A) < \/2. Thus, BHAL is valid for a certain 
range of 7 which is same as what has been reported in [22] i.e. 0.159 < 7 < 0.225. Hence, for this particular 
range of the Barbero-Immirzi parameter (7), we have the microcanonical entropy of an IH given by 

3 

Saw = S B h - ^ lo S s bh H (29) 

where Sbh = Ad/4£* is the Bekenstein- Hawking entropy and the smaller terms are neglected. From 
cq.(29), one can clearly observe that the universal coefficient in front of the logarithmic correction is —3/2 
as usual[23, 28]. But, from eq.(25), it is revealed that while expressed in terms of k and N, there is a 
fundamental effect due to the presence of the topological defects, the punctures, playing the role of a new 
'quantum hair' first proposed in [19]. Also we have an appropriate bound on the permissible values of 7, 
considering the BHAL as an external input as pointed out in [22]. 



V. DISCUSSION 



As mentioned in the Introduction, the LQG formulation of a QIH in terms of a CS theory coupled to spins 
directly involves the use of the quantum parameters ('hairs') k and N. Strictly speaking, if the classical area 
A c i is taken to be a hair, the Barbero-Immirzi parameter 7 (the coefficient of a topological contribution to 
the classical action from the Nieh-Yan invariant [31]) is considered as an independent coupling parameter, 
then the only new parameter that has appeared in the quantum theory is the number of punctures N. 
Equivalently, one can, as we have in this paper, take k and N to be the parameters characterising the 
quantum theory. It is significant that the requirement that the microcanonical entropy of a QIH yields the 
area law for large k , N immediately restricts the Barbero-Immirzi parameter 7 to lie in a specific interval 
on the real line, while also yielding the subleading logarithmic correction derived in earlier literature with 
the universal coefficient —3/2. Thus, the complete characterization of a QIH is given by two independent 
parameters, namely, k and N. 

The departure from recent literature is the direct link established in this paper of these results and inter- 
pretation within the CS formulation of a QIH. The point here is that the computation of the microcanonical 
entropy of a QIH is not merely a combinatoric exercise in statistical mechanics of an ideal gas of punctures 
with some restrictions gleaned out of LQG as in some part of the recent literature [19, 32-37], with at best a 
weak link to the complete formulation as a CS theory that has been available for years. Rather, the original 
formulation yields an understanding that is complete as a quantum theory. 

An extension of the foregoing analysis to Trapping horizons and Entanglement entropy of matter and 
radiation fields in their vicinity, would amount to truly new physics of radiant quantum horizons beyond 
general relativity. This has within it the potential to surpass, because of its firmer quantum geometric 
underpinning, recently proposed speculative ideas of a somewhat ad hoc nature based on semiclassical 
analysis [38, 39], about how classical horizon geometry must change (become a 'fuzzball' or a 'firewall') 
so as to allow the existence of well-defined scattering amplitudes for quantum matter fields. Whether 
or not such a structure emerges from LQG in an appropriate limit is not known at the moment, since 
there is still no complete quantum geometric analysis of Hawking radiation and its various conundra from 
an LQG standpoint. The problem at hand involves a quantization of the kinematical phase space of a 
Trapping/Dynamical horizon along the lines of [1, 2] which may not be technically so simple as a QIH 
because of the transient behaviour inherent in the geometry. However, perhaps a perturbation of the CS 
description by a set of appropriately chosen matter field operators might serve as a first approximation to 
the problem. The strength of the LQG approach lies in the transparent manner in which the CS symplcctic 
structure emerges from the boundary conditions in the incipient formulation [2], without having to rely on 
conjectured results. The relation of the CS Hilbert space to the conformal blocks of the WZW model on 
the QIH is also not a matter of conjecture for large black holes. Thus, the LQG analysis of a QIH geometry 
throws up holographic structures as emergent, without any prior notion that they have to be there. One 
expects that generalizations to Trapping horizons to centre around this theme so as to reap the benefits of the 
4 dimensional gravity - 2 dimensional conformal field theory link found in ref. [20]. Note however that the 
thermal stability of radiant trapping horizons which approach an equilibrium QIH has been discussed within 
the LQG framework yielding a criterion of stability involving the equilibrium mass and the microcanonical 
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entropy of the QIH in ref. [40 42] . 
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